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Abstract 



Possible sources of uncertainties in the calculations of the internal conversion 
^ . coefficients are studied. The uncertainties induced by them are estimated. 
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a) '. 1 Introduction 

The internal conversion coefficients (ICC) serve many years as an indispensable tool to as- 
sign spin-parity quantum numbers of excited nuclear levels. There exist widely distributed 
tables [I], |2|, 0] of ICC covering broad region of transition energies, atomic numbers from 
Z=30 H, |] or 10 @ up to Z=103 or Z=10A [| § and transition multipolarities Ml 
- M4, El - E4. Very recently, tables of ICC for the superheavy elements (104<Z<126) 
appeared [§J. 

When computing ICC, however, various approximations are applied. Some of them 
are of purely theoretical origin, e.g. the calculations are done only for the first non- 
vanishing order of the perturbation theory. Other ones are connected with improper 
description of physical reality, e.g. the ICC is calculated for a free atom whereas the real 
source atom is usually chemically bound in some compound or embedded into a solid. 
In this work we systematize these particular sources of uncertainties and estimate their 
effect on the accuracy of the calculated ICC. 

*Main part of this paper was presented at the meeting of Decay Data Evaluation Project, Braun- 
schweig, May 2000 

^e-mail: rysavy@ujf.cas.cz 



1 



2 Basic formulae 



The basic formulae for the ICC together with its derivation may be found in many works. 
We present here those from the monography ||. The ICC for the transition of the 
multipolarity tL and the atomic shell i is given by 

Here 



n( ML) = / ^(0) (r)M ^ r)rfr; (2) 

Jo 

n[ EL ^ = / [LF K (r)h L (ur) + ^(r^x-iM^r. (3) 
Jo 



and 



^Kfi r ) = (^i + K f )(u Ki v Kf + u Kf v Ki ), (4) 
F K} (r) = (u Ki u Kf + v Ki v Kf ), (5) 
^ 1) ( r ) = ( K i ~ K f ~ L)v Ki u Kf + (Ki -K f + L)u Ki v Kf . (6) 

^j'umi,j2,m 2 ^ S Clebsch-Gordan coefficient, a is the fine structure constant, h t (x) is the 
spherical Hankel function, uj is transition energy, u K and v K are small and great, respec- 
tively, component of the solution of Dirac equation for electron. 
For the transition of magnetic multipolarity, the quantity 

T(ML) = I d 3 R Jnud(R) L g^(R) Yju 
Kf J d 3 R J nucl (R) t j L (uR) Ytu 

where 

gf{R) = h L (uR) [ dr j L {ur) F K {r) - j L {uR) f dr h L {ur)F K {r) (8) 
Jo Jo 

is the only one which containes nuclear variables. It describes "intranuclear conversion" 
(see Sect. |3.7| ). For the electric transitions, the analogical quantity is given by a 

more complicated formula but it has similar structure. 

From the formulae above we see that the ICC depends on four 'parameters': tran- 
sition energy u, transition multipolarity tL, initial electron state (i.e. atomic shell), and 
atomic number Z. For the atomic number, the dependence is purely implicit - via atomic 
potential which then forms the electron wave functions u K and v K . The other quantities 
enter into the formulae explicitly but, moreover, they also affect the result indirectly. 
E.g. the transition multipolarity determines the allowed final electron states due to the 
conservation of angular momentum. It is obvious that to draw some information on the 
shape of the dependencies directly from the above formulae is almost impossible. The 
only way is to evaluate ICC for various parameters and find out the dependencies from 
the numerical results. Analogically, the effect of various approximations and/or model 
imprefections may also be studied via numerical calculations. 
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3 Sources of the ICC uncertainties 



When we put a question about the precision of the theoretical ICC we must first specify 
the possible sources of their uncertainties. These may be very roughly divided into two 
classes — those stemming from theory and those given by our insufficient knowledge of 
some quantities entering the calculations. Into the first category there can be included 
atomic model, nuclear deformation and nonsphericity in general, and higher order effects. 
The second one consists of atomic binding energies, isotonic effect, and chemical effects. 
Somewhere in-between is the question of correct accounting of nuclear structure, i.e. the 
intranuclear conversion. And, finally, a serious error may be introduced by incorrect 
interpolation in existing tables, especially in case of total ICC. 

In the following, we address ourselves to the particular items listed above. We 
summarize all the known facts and try to give a realistic estimate of the uncertainties 
induced. 



3.1 Atomic model 

As seen from Eqs. (§ - H), the electron wave functions (u K and v K ) are needed to calculate 
an ICC. The evaluation of these wave functions is based on some suitable model of atom. 

When we skip the 'ancient' calculations with the model of point nucleus without 
screening, the first more realistic ICC values for the K and L shells were obtained 0-0 
using the Thomas-Fermi-Dirac atomic model. Later on this model was abandoned in 
favour to physically more founded model of Hartree and Fock with the Slater exchange 
term (HFS). This atomic model was employed to provide the first tables of the N-shell 
ICC ||. A little sooner in Ref. a modified HFS method where the Slater exchange 
term is weighted by a factor of 2/3 (so called HF|S method) was utilized to get the ICC 
for the shells K, L, and M. 

In the late 70', the most extensive tables |2| appeared where the pure HFS model 
was applied. Approximately at the same time the first ICC based on the true Hartree-Fock 
(HF) model (without the Slater exchange term) were calculated [|K], [11]. To complicate 



the things yet more, it is not clear whether the 'hole' in the atomic shell remaining after 
the conversion electron should be taken into account or not. 

Evidently one should hesitate which ICC are the "real" ones and what error can be 
made by using other ones. Band et al. [0 performed extensive calculations of ICC in 
both models HFS and HF|S and for both variants 'hole' and 'no hole'. They evaluated 
large sets of ICC for Z=30, 60, 90, conversion electron energy between 10 and 400 keV, 
multipolarities Ml - M4 and El - E4, and K, Li, L2 and L3 subshells. Then they 
constructed the quantities which may characterize the uncertainty caused by the atomic 
model, 

r ICC(hole) 

Ai = jrr( 7 ' - 1 x 100, 9 
ICC [no hole) 

JCC(HF\S) , , , 

= ' ICC ( HFS) ~ 11 X 100 ' (10) 
where the symbol in parentheses indicates the model used. These quantities vary with 
multipolarity - in general, their absolute values increase with L. They increase with 
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decreasing energy. For some combinations of subshell and multipolarity they do not 
exceed 0.5%, for the remaining ones there are graphs in [TL2 



In [pjj] , an extensive comparison of HF modelQ ICC with experimental ones was 



performed. Altogether 57 experimental data items for 8 transitions in six different isotopes 
were examined. A sensitive statistical test revealed better agreement of the theory with 
experiment for the HF ICC than for the HFS ones. The differences between the two sets 
of theoretical ICC varied from <1 % on the K-shell, 2-3 % on the L-subshells up to ~30 % 
on the outermost shell. 

More recently, there appeared [0 the 'true' HF model ICC, i.e. those where even 
the free-electron wave functions were evaluated without the approximation of the exchange 
term. 15 transitions in various isotopes were studied and better agreement of the theory 
and experiment was confirmed. 

Tables of ICC in the HF model are not available at present. However, as stated 
above, programs exist and it is possible to calculate the HF (i.e. physically best) ICC for 
particular cases. As for their uncertainty, only that corresponding to numerical accuracy 
(mostly <0.5 %) may be expected. For the other atomic models, it is difficult to give 
some general error estimate. It might be the difference from the HF ICC. Anyway these 



quantities are known for several particular transitions only (e.g. [IT], And in any 

case, such an 'error' would be a systematic rather than statistical one. 

Recommendation : For analysis of precise experimental data, use the HF model based 

ICC. 



3.2 Nuclear deformation and nonsphericity 

The formulae for conversion coefficients are based on the assumption that the atomic field 
is spherically symmetrical. 

It has been proved long time ago |15| that the closed shells form a symetrical atomic 
field. With several exceptions ( ||16|| , p. 74), there is only one unclosed shell in every free 
atom. This is the valence shell which, in most cases, is the outermost one. (Even if there 
are two unclosed shells, they always belong to the outer ones.) The calculated ICC is in 
fact the mean over all projections of the angular momentum. This evidently brings an 
error in the case of unclosed shells. If the actual occupation number for the particular 
shell, i, is n act and the maximum possible one is n Q , the conversion coefficient should be 
rescaled as 

aAact) = — x aAcalc), (11) 
n 

where a;(ca/c) is the calculated ICC. It greatly reduces the error caused by the absence of 
one or several electrons. The error due to the asymmetry remains but it can be negtected. 
The asymmetry is in fact a slight perturbation which affects almost exclusively only the 
particular unclosed shell. As stated above, it concerns the outer shells and there is a 
greater uncertainty due to the chemical environment (Sect. |3.6| ). 

1 In fact, only the bound electron wavefunctions are calculated in HF model there. The wavefunctions 
of the emitted electrons are generated as previously is HFS model. It was proved that the free electron 
wavefunctions are much less sensitive to the model details (e.g. |i~3|]). 
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Another possible source of asymmetry can be a deformed atomic nucleus. The 
effect of nuclear deformation on the ICC was studied in [17, 18, [BJ. The authors [17, |T/g 



assumed a nucleus with high quadrupole deformation and used perturbation series to 
compute contributions to ICC. Such calculations were done for ^Dy, 1 ggEr, ^Yb, and 
1 ygOs nuclei exhibiting high quadrupole moments. The comparison of the results with the 



ICC evaluated in a usual way 

-2 



40" 



70. 



In 19 



i.e. for spherical nucleus) revealed a negligible change of 
also the effect of magnetic dipole moment was investigated and found 



to be negligible. 

It should be kept in mind that the transitions in deformed nuclei are often hindered. 
In such cases it may be necessary to account for the intranuclear conversion (Sect. |3.7| ). 

Recommendation : If using the conversion coefficient for a subshell which is not 
completely filled do not forget to rescale it following Eq.(Pl). If studying a transition in 



a deformed nucleus, think on possibility of intranuclear conversion. 



3.3 Higher order effects 

At present, all ICC are evaluated up to the first non- vanishing order of perturbation 
theory, i.e. the first order for the gamma emission and second one for the emision of 
conversion electron. The Feynman diagrams corresponding to the next nearest orders 
are presented in Fig.Q. The contributions of individual graphs to ICC were estimated by 
several authors, e.g. p0f-|30|,H. 





Figure 1: Feynman diagrams for internal conversion up to the fourth order and for 7-ray 
emision up to the third order of perturbation theory. 

Vacuum polarization. It is the only higher order effect which can be taken into 



acount precisely by a change of atomic potential. It was studied in [2£|. The authors 
found small influence on ICC, namely a change of <0.1 % for Z=26 and <1 % for Z=80. 

Vertex correction. Contribution of vertex correction to ICC was never investigated. 
Pauli et al. |25| stated that Hager and Seltzer |24| studied it. However, the graphs 
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presented in |24| correspond to 'electron bridge' - graphs (e) and (f) in Fig.|T]. 



Self-energy. Up to our knowledge the possible contribution of this term to the ICC 
was never examined. 

Bridges. These corrections were extensively investigated by Krutov et al. ||20||— 123| . 



Unfortunately, there are almost none numerical estimates. The up to 20 % effect [22 
of the graph (d) in one specified case seems unrealistic (during the calculations, a lot of 
terms which can contribute were neglected). Hager and Seltzer pi] stated about 5 



efect of (e) and (f) on the L-subshell ICC in ^W, 1 ggTm, and ^Dy. Also here, several 
important terms were neglected and, moreover, the Pauli principle was not fulfilled in 
their calculations. On the other hand, Mayol et al. [^7] found the electronic bridge effect 
of 0.0047 % in 77 Ir atom which is stated to be the maximum one from 11 isotopes studied. 
Hinneburg et al. [J2l| reported extremely small contribution of 10~ 9 . In the later works 
|29|, [30|[ they concentrate to the transition with the energy of 77 eV in 235 ™U where the 



92 

graph (f ) could raise the photon decay constant by a factor of ~10 5 and thus suppress the 
ICC. Their calculations, however, cannot be verified since the ICC of such a low energy 
transition is unmeasurable at present. Band et al. |bl admit the effect of 'electron bridge' 



for very low transition energies but do not present any concrete estimates. 

Anyway, the above approach - the evaluation of contributions of individual graphs 
- is not correct. The right method is first to derive the matrix elements corresponding to 
particular graphs and add them together. Only then the modified integrals (0,0,0) can be 
applied in the final formula ([!]). This is such a difficult task that nobody even dared to 
perform it up to now. 

Recommendation : For the transition energies not too close to conversion threshold, 
assume the uncertainty caused by neglection of higher orders to be about 1 %. This 
estimate follows from analogy with electron Coulomb scattering using the formulae in 



[ 3j| . As for the very low transition energies, we can say nothing at present. 



3.4 Atomic binding energies 

The binding energies, E^, of the atomic electrons enter into the ICC via the relation 

Ekin,i = E tr — Eb t i, (12) 

where Eki n ,i and E tr are kinetic energy of conversion electron and transition energy, re- 
spectively. In Eq . ([12|) , the recoil energy as well as work function of the source material 
are neglected. The Eki n ,i specifies the conversion electron wave functions appearing in 
formulae (|]-|6|). 

The binding energies are presented in [B2[ with the typical error of ±1 eV, some- 



times 2 or 3 eV. (An exception are several transuranic elements where the stated error 
amounts in some cases up to ±18 eV.) To estimate the uncertainty in ICC caused by those 
uncertainties in E^ we performed calculations with various values of binding energies. In 
particular, we evaluated ICC for 26Fe, 55CS and 87Fr and the subshells (K,M 3 ), (K,N 5 ) 
and (K,Os), respectively. This setting was chosen to cover low, medium and high Z and 
always one inner and one outer shell. The transition energies of 30, 50, 100, 200, 300 and 
500 keV were used where applicable. Finally, the ICC were evaluated always for three 
binding energies, E -10 eV, E , and E +10 eV where E Q is the binding energy f32fl . (Note 
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that the binding energies - including those for the inner shells - may also be influenced 
by the chemical state of the converting atom. These so called chemical shifts, however, 
are rather small and our interval of 20 eV covers them all right.) 

Table 1: Internal conversion coefficients for 55CS, shell K, transition energy of 50 keV and 
three binding energies. 



mult. 


Eb-10 eV 


E a 


E +10 eV 


A b 


El 


1.2320(0) c 


1.2318(0) 


1.2317(0) 


.02 


E2 


7.0612(0) 


7.0563(0) 


7.0513(0) 


.07 


E3 


2.6671(1) 


2.6632(1) 


2.6593(1) 


.29 


E4 


9.2835(1) 


9.2622(1) 


9.2408(1) 


.46 


Ml 


5.8912(0) 


5.8912(0) 


5.8911(0) 


<.01 


M2 


9.4640(1) 


9.4617(1) 


9.4594(1) 


.05 


M3 


6.1914(2) 


6.1863(2) 


6.1811(2) 


.17 


M4 


3.1418(3) 


3.1369(3) 


3.1320(3) 


.31 



E =36189.9 eV 



b Maximum change in percent, A = (2|£ _ 1) x 100 
c x(y) means xxl0 y 

As a measure of the uncertainty we take the quantity A = (222s _ 1) x 100 where 
max and min are the maximum and minimum, respectively, ICC from the three ones 
corresponding to different binding energies. For the innermost (K) shell, the results are 
following. In almost all cases, A is of order of 10~ 2 of percent, up to 0.1 %. The exception 
is the 50 keV transition in 55CS (see Tab.^). For high multipolarities, the changes there 
reach almost 0.5 %. The reason is simple - we are closer to threshold than in the other 
cases. (When approaching the threshold, the kinetic energy decreases and its change plays 
more important role.) For the outer shells, the changes are even less. 

The change of binding energies used here (20 eV) is relatively large. If the ex- 
perimental error of binding energy is less (which is usually the case), the corresponding 
uncertainty in ICC is also less. 

Recommendation : If the binding energy is known with a precision better than say 
±5 eV and the transition energy is high enough (20 keV or more above the threshold), 
the uncertainties caused by it in ICC are negligible. This is true even when we keep in 
mind the possible chemical shifts of binding energies. For the transition energy very close 
to threshold, no general rule can be given and the problem must be solved by calculating 
ICC for different binding energies. 

3.5 Isotopic effect 

The mass number, A, does not enter the formulae for ICC explicitly. Nevertheless it 
determines the nuclear radius 

rnuci = r x A 1/3 . (13) 
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This implies that A affects the atomic potential which determines the electron wave 
functions. We describe the atomic nucleus by the Fermi distribution of charge 

<*■> - TT^W-y (14) 

where c and r nud are tabulated in [[33|] and t is found from the condition p{r nud ) = 0.05p - 
We assume that r nud in corresponds to A of the most abundant isotope for 



every element. To obtain ICC accordant with another mass number, A', we rescaled the 
nuclear radius by the ratio (A' /A) 1 ^ 3 . The calculations were performed for the same set 
of elements, transition energies, shells and multipolarities as that in Sect.( j3~4[ ). ICC were 
evaluated for three different mass numbers, A and A ± 10, in 55CS (,4=133) and 87R: 
(,4=223) and for two mass numbers (57 and 67) in the case of 26F e - 

The results are following: In iron and cesium, only very small changes (up to ~0.2 %) 
were observed for the K-shell ICC. For the outer subshells considered, the changes were 
significantly lower. In the francium, the changes turned out to be higher (see Tab.^). 
Generally, the changes do not depend to transition energy. This indicates that they are 



Table 2: Internal conversion coefficients for 87Fr, shell K, transition energy of 200 keV 
and three mass numbers. 



mult. 


,4=213 


,4=223 


,4=233 


A a 


El 


7.0255(-2) 


a 7.0251(-2) 


7.0250(-2) 


0.01 


E2 


1.6199(-1) 


1.6197(-1) 


1.6196(-1) 


0.02 


E3 


3.8503(-l) 


3.8490(-l) 


3.8473(-l) 


0.08 


E4 


9.6301(-1) 


9.6236(-l) 


9.6176(-1) 


0.13 


Ml 


1.7636(0) 


1.7586(0) 


1.7538(0) 


0.56 


M2 


6.6561(0) 


6.6465(0) 


6.6379(0) 


0.27 


M3 


1.7317(1) 


1.7295(1) 


1.7277(1) 


0.23 


M4 


4.2730(1) 


4.2689(1) 


4.2642(1) 


0.21 



See comments on Tab.|I] 



caused mostly by slight alterations of the bound-electron wave functions. As could be 
expected, the K-shell is affected most of all. Moreover, the magnetic transitions are more 
sensitive than the electric ones. 

Recommendation : In most cases, the uncertainties caused by isotopic effect are 
negligible. In case that data on some extremely heavy or extremely light isotop of a 
given element are studied, allow for uncertainty of up to 0.5 % to be on the safe side. A 
better alternative is to evaluate ICC directly for the given A which eliminates this kind 
of uncertainties completely. 



3.6 Chemical effects 

The conversion coefficients are always evaluated under the assumption that the source 
atom is an isolated free neutral atom. In practice, the experimental ICC are obtained from 
a radioactive source where the radiating atoms are bound in some structure - molecule 
or crystal. 
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There were attempts to account for the source structure. The Wigner-Seitz bound- 
ary condition introducing solid state effects was applied in calculations |[27|| . The ICC 
were evaluated for 11 transitions in various isotopes. It turned out that the K-shell and 
total ICC were almost the same as those evaluated in the free atom approximation. (As 
for the other shells, no information is given in [27].) In [34], the ICC for the free atom 



were compared with those calculated using the Xa-SW description of a Tc compound. 
The difference between the two total ICC was of order of 10 _1 %. In [[35], the Xa-SW 
calculations for some tellurium compounds gave electron densities near the Te nucleus. 



The free-atom ICC were then scaled as lCC(real) = \ip Xa ~ MS [r nud )\* / \ip J 1 ^{r nud )\*x 
ICC(/ree), where |^| 2 's are the bound electron densities at the nuclear edge. Unfortu- 
nately this attepmt failed to bring the theoretical ICC to an agreement with experimental 
ones for several Te compounds. 

Another possibility is to approximate the 'real' ICC by their values evaluated for 
free ions. The atom bound e.g. in a molecule certainly is not a free ion. In spite of this the 

2 gfBi, and 1 |gPr, respectively) 



ICC for free ions were calculated 36 - 39 



(for iTe, 



83 K 



and it turned out that the experimental outer-shell ICC ratios lay within the range of 
theoretical values for reasonably ionized atoms. This gives us an argument to use the 
dispersion of ICC calculated for several ionic states as a measure of uncertainty caused 
by the chemical bond. 

First, we performed extensive calculation of ICC in 1 f|Cs over a broad region of 



ionizations (from -1 up to +9) 
shells, and the energies of 30, 
calculation is described in [FT? 



The multipolarities Ml - M4 and El - E4, all atomic 
100, 200, and 500 keV were treated. The method of 



3.6 



3.2 



2.8 



X 
c. > 



2.4 - 



2.0 - 



1.6 



-2 



Shell: 

Energy: 

Multipolarity: 



30 keV 

El 



evaluated 

corrected for partial occupancy 



- 1 



ionicity 



Figure 2: ICC on the Oi-subshell of 55G3 in dependence of ionization of the atom. With 
the exception of the ionization +8, the Oi-subshell is fully occupied. 



As a measure of uncertainty we used again the quantity A introduced in Sect.|3^4 
(Tab.Q) but now the maximum and minimum is always taken over all ionizations. The 
magnitudes of A were as follows: <0.2 % (for the shells K, Li-Lg, M1-M3), ~0.3 % 
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(M4, M 5 ), ~0.6 % (Ni), increasing up to 10 % (for N 5 ), and several tens of percent 
(O-subshells). For the total ICC, A was mostly ^0.3 %, exceptionally ~0.5 % (for Ml 
transition, energy 30 keV). The results are not too dependent on both energy and mul- 
tipolarity. For illustration there are results for the subshell M2 in TabfJ. The changes of 
ICC on the Oi shell shows Fig.Q. 

Table 3: Uncertainties (in percent) of ICC on the M 2 -subshell of 1 || Cs caused by chemical 
binding. 





E 7 [keV] 






mult. 


30 


100 


200 


500 


El 


0.10 


0.09 


0.09 


0.07 


E2 


0.07 


0.07 


0.11 


0.08 


E3 


0.07 


0.07 


0.07 


0.11 


E4 


0.17 


0.08 


0.06 


0.07 


Ml 


0.10 


0.13 


0.10 


0.08 


M2 


0.08 


0.07 


0.10 


0.07 


M3 


0.20 


0.09 


0.11 


0.06 


M4 


0.22 


0.10 


0.07 


0.07 



To get more information on the 'chemical' uncertainties, we extended our study to 
more elements. In particular 87 Fr - together with 55CS - belongs to alcaline metals, i.e. 
those having one electron above the noble-gas core. Then 3sBr, 53 I, and §sAt - here, 
one electron lacks to noble-gas-like closed shells. And finaly 26Fe, 62SH1, and 76 0s - the 
transition metals. This set covers various configurations of the valence shell. Moreover, we 
added the transition energies of 50 and 300 keV to suppress possibility of some unobserved 
changes of the uncertainty. 

The results confirmed the conclusions made from the case of 55 Cs. The changes are 
of tenths of percent for the iner shells and increase to tens of percent for the outermost 
shell. The changes (uncertainties) of the total ICC are presented in TabfJ together with 
the ranges of ionization used for particular elements. We see that the uncertainty of total 

Table 4: Maximum changes of total ICC over all the ionizations given 



ionization change 





from 


to 


[%] 


26 Fe 





+4 


0.7 


35 Br 


-1 


+4 


0.5 


53I 


-1 


+4 


0.2 


55 Cs 


-1 


+9 


0.5 


62 Sm 





+5 


0.1 


76 0s 


-1 


+5 


0.1 


85 At 


-1 


+4 


0.1 


87 Fr 





+8 


0.2 



ICC due to the chemical effects decreases with increasing atomic number, Z. This is not 
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surprising - the higher Z, the more atomic subshells. And only the lower ones (which are 
less affected) contribute significantly into the total. 

Recommendation : If no information on the structure of the source exists, the quan- 
tities presented in Tab.f| may be used as very conservative estimates of the error. On the 
other hand if some information is available, the present values may be modified (usually 
decreased) or neglected. 



3.7 Intranuclear conversion 

The intranuclear conversion (also called nuclear structure effect, penetration effect) is 
responsible for that part of ICC which corresponds to the integral (0). Physically it 
is caused by non-zero probability that the electron penetrates into the nucleus and the 
conversion process takes part there. 

In principle, the nuclear currents J nU ci(R) in Eq.(0) may be determined from a 
suitable nuclear model and the integral then may be evaluated. Such an ICC, however, 
would be model-dependent. Moreover, the reliability of nuclear models is not sufficient. 
Therefore, two approximations were introduced. The first one, the Rose's 'no penetration' 
model P, pED[], completely neglects the contribution (|7|). The second one, the Sliv's 'surface 



current' model |Q assumes that the nuclear transition currents are localized at the nuclear 
surface which enables the evaluation of the integral ([7|). 

Note that the both approaches are approximations only (even when each of them is 
supported by some physical basement). A convenient way how to manage the contribution 



of the intranuclear conversion was developed by Pauli ||^] . He expressed the ICC involving 
the contribution (0) as 

aa(ML) = a?(ML) x (1 + b u \ + b 2l \ 2 ) (15) 
ai(EL) = a>i(EL) x (1 + a u r] + a 2i r] 2 + a 3i r]( + a 4i £ + a 5i £ 2 ) (16) 

where a° is the ICC on the subshell % calculated in the 'no penetration' model, the formulae 
for a's and b's contain only electronic variables. All nuclear variables from Eq.(0) are 
constricted into one (A) and two (£, 77) parameters for the ML and EL multipolarities, 
respectively, which, in turn, contain no electronic variables. These parameters should be 
determined by comparison of the theoretical and experimental ICC using e.g. the program 
or, in principle, they may be calculated in frame of a specific nuclear model. 
The possible effect, i.e. the magnitude of the polynomials in parentheses in Eqs. 
I^,r^), was extensively studied in 0, [0J. For non-hindered transitions, the effect - the 



change of a° - varied from negligible ~10 -2 % (for low Z and all multipolarities except 
Ml and M2) to 10 % (for high Z and multipolarity Ml). For the hindered transitions, 
the situation is different. In ||, a lot of various transitions with the hindrance factors 
F w (related to the Weisskopf single particle estimate) from 10 up to 10 16 were studied. 
It turned out, that the effect fluctuates between a few percent change up to a change 
by a factor of 20. Unfortunately, there is no universal connection between the hindrance 
factor and the magnitude of the effect. As an example, the 58 keV El transition in ^Hf 
has F\y = 10 16 and the effect reaches the factor 2.5 |45| |. On the other hand, the 84 keV 



El transition in 2 f*Pa has Fw = 9xl0 5 only and the ICC changes up to 20 times |46[ 
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(Note that the change due to the intranuclear conversion is different for different atomic 
subshells. When present, the effect changes the ICC not only for the inner shells but also 
for the outer ones ||7|, |48|.) 

It is necessary to emphasize here that the intranuclear conversion itself brings no 
uncertainty into the theoretical ICC. The theoretical ICC are calculated using one of the 



above models || |4ip . In many cases - when the effect of nuclear structure is small - 
the theoretical ICC enabled the determination of nuclear characteristics (spins, parities) 
fairly well. However, when treating precise experimental data or hindered transitions, 
the intranuclear conversion must be taken into account. This means that the ICC, cti, in 
the form of Eqs.(Jdj, Tb\) must be used instead of a® and the parameters A, rj, £ should be 
fitted [^3J into the experimental data. The parameters are determined together with their 
errors (standard deviations) but these errors are result of the fit and they do not affect 
the theoretical ICC 

Recommendation : If there is a discrepancy between the experimental and theoretical 
ICC and/or ICC ratios which cannot be explained by higher multipolarity admixture, 
always expect the intranuclear conversion and do fit the parameters A or £ and rj. 



4 Danger of interpolation 

If there is no possibility to calculate ICC directly for the transition in question, one 
is forced to use ICC tables. The tables are also usually used for preliminary estimate 
of effects which can be expected in the detailed analysis. The tabulated ICC are always 
presented in some more or less dense energy mesh. When studying a particular transition, 
an interpolation is necessary. 

It is generally believed that the dependence of ICC on the transition energy is close 
to linear one in the log-log scale (i.e. log ICC ~ log-E^). If so, the interpolation would 
make no problems. Unfortunately, in several cases the situation is not so simple (see e.g. 
graphs in ^S|). For low-energy transitions a very non- monotonous dependence of ICC on 
energy was found in [|50|j . 

For the total ICC, the interpolation is strictly forbidden in the vicinity of the thresh- 
old of some atomic shell. As demonstrated in Fig.|3] there is a discontinuity in each point 
E tr = Ebj where Ej,^ is the binding energy of some atomic subshell. The magnitude of 
this discontinuity depends on the transition multipolarity and on the particular subshell. 
(E.g. for Ml, the discontinuity is greater for S1/2 (i.e. K, Li, Mi etc.) subshells and less 
for pi/2 and p%/2 subshells. For E2, it is just opposite.) Applying an interpolation in such 
a case would fetch in a very serious error. It is a pity that this error was not avoided in 
pifl (in graphs of total ICC, appendix F). 
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Figure 3: Depencence of total ICC on transition energy for the E2 transition in element 
Z=126. (Note that no discontinuity is observed at the si/2-shell thresholds since the 
contribution of these shells into the total ICC is very small in this case.) 
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5 Conclusions 



We have estimated the uncertainty of the theoretical ICC relevant to their most frequent 
application, i.e. the spin-parity determination of excited nuclear states. Until the un- 
certainties of the experimental ICC and/or their ratios are better than 10 %, any of the 
three tables ]I], ^ |3| can be utilized for the comparison with the experiment and their 
uncertainties may be neglected. 

When the accuracy of the experimental data is of a few percent, the uncertainty 
of the theoretical ICC as examined in this work should be included into the analysis 
which, in turn, leads to more realistic error of the final results (e.g. the multipole mixing 
parameters). As for the adding the uncertainties, we recommend to apply the quadratic 
addition as in the case of the statistical errors. We admit that the mentioned uncertainties 
are not of the pure statistical origin but the same argument as for the statistical errors 
can be used - it is improbable that all deviations assemble into the same direction. 

To circumvent the question of errors caused by the atomic model which is not clear 
enough, the ICC for the cases of precise experimental data shoud be calculated directly 
for the particular transition energy using the HF atomic model (see Sect {3 . 1| ) . 

Keeping in mind the possibility of nuclear structure effect (Sect j3.7|) in hindered 
transitions, more than one experimental quantity (ICC and/or ICC ratio) shoud be fetched 
into analysis whenever possible. 

The interpolation in the tables of total ICC should be avoided. This is due to 
discontinuities of total ICC at the points of particular subshell's threshold. For lower 
transition energies, even the dependence of subshell ICC on the energy becomes non- 
monotonous and the interpolation must be done with a great care. 
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